Fermi edge singularity in a non-equilibrium system 
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We report exact results for the Fermi Edge Singularity in the absorption spectrum of an out-of- 
equilibrium tunnel junction. We consider two metals with chemical potential difference V separated 
by a tunneling barrier containing a defect, which exists in one of two states. When it is in its 
excited state, tunneling through the otherwise impermeable barrier is possible. We find that the 
lineshape not only depends on the total scattering phase shift as in the equilibrium case but also on 
the difference in the phase of the refiection amplitudes on the two sides of the barrier. The out-of- 
equilibrium spectrum extends below the original threshold as energy can be provided by the power 
source driving current across the barrier. Our results have a surprisingly simple interpretation in 
terms of known results for the equilibrium case but with (in general complex-valued) combinations 
of elements of the scattering matrix replacing the equilibrium phase shifts. 
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Developments in the fabrication and manipulation 
of mesoscopic systems have allowed detailed and well- 
characterized transport measurements for a large range of 
devices including quantum pumps, tunnel junctions and 
carbon nanotubes. It is often the case that such measure- 
ments explore non-equihbrium effects particularly when 
the potential difference is dropped across a narrow po- 
tential barrier or over a short distance inside the metallic 
region P, 0, Q. While there is often a very good the- 
oretical description of much that has been observed for 
systems close to equilibrium, the theoretical picture for 
systems out of equilibrium is less clear with fewer estab- 
lished theoretical results. 

A natural point to start, when looking for a descrip- 
tion of non-equilibrium effects in many-electron systems 
is the Fermi Edge Singularity (FES), which is one of the 
simplest non-trivial many-body effects. The FES is char- 
acteristic of the response of a Fermi gas to a rapid switch- 
ing process. Initially it was associated with the shape of 
the absorption edge and spectral line found when a core 
hole is created However, it turns out to be a generic 
feature of a Fermi system's response to any fast switch- 
ing process and reflects the large number of low-energy 
(particle- hole) excitations which exist in Fermi liquids. It 
has also been shown to be related to Anderson's orthog- 
onality catastrophe |^ and can be used to reformulate 
the Kondo problem in terms of a succession of spin flips 
which are treated as the switching of a one-body poten- 
tial between two different values Q. 

We consider a system at zero temperature with two 
Fermi surfaces separated by a barrier with a potential 
difference (bias) V applied across the barrier (see Fig. 

The barrier contains a defect, which exists in one of 
two states with energy separation Eq. Tunneling through 
the barrier is assumed to be possible only when the de- 
fect is in its excited state. We compute the absorption 
spectrum close to the threshold at = Eq — 'Rc{IS.{V)) , 
for frequencies (tj — tJo) ^ foj where is of order the 




FIG. 1: Energy levels in an idealized device to demonstrate 
the out-of-equilibrium FES. The scattering potential for elec- 
trons is characterized via the 2x2 matrix, S{e, t), connecting 
scattering states in the two wires for particles with energy e. 
S — OT 5"^ depending on whether the defect is in its ground 
(,g) or excited (e) state (with excitation energy Eo). is the 
identity matrix and is an arbitrary unitary matrix. Sfi and 
correspond to the reflection and transmission amplitudes 
respectively. 



bandwidth and Re(A(l/)) is real part of the combined 
energy shift of the two Fermi seas when the defect is in 
its excited state. (A(y) is complex for non-zero V on 
account of the dissipation in the system.) 

Using an approach based on that of Nozieres and de 
Dominicis (ND) 0, we solve exactly for the asymp- 
totic behavior of the absorption spectrum in two limiting 
cases: (uj — luq) ^ V and {uj — loq) <C V. Our results have 
a simple interpretation in terms of generalized (complex) 
phase-shifts at the Fermi energy. Typical lineshapes for 
the case {uj — ujo) ^ V illustrating the dependence on the 
reflection amplitudes and phases are shown in Figure [3 

Our treatment of the problem is based on that of 
Muzykantskii and Adamov (MA) for the statistics of 
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charge transfer in quantum pumps, which uses the re- 
lation between the many-particle response to the chang- 
ing one-body potential and the solution of an associated 
matrix Riemann-Hilbert (RH) problem This prob- 
lem was also addressed perturbatively and using the ND 
approach in 0, 0| , although the results in led the au- 
thors to question the validity of the ND approach of ^3 
(sec also Our solution shows clearly that the ND 

approach is valid, with the earlier difficulties probably 
associated with an incomplete analysis of the matrix RH 
problem associated with their singular integral equation. 

We characterize the scattering at the interface between 
the two subsystems via the unitary 2x2 matrix, S{e, t), 
connecting scattering states in the two wires for parti- 
cles with energy e. This takes one of two values S^, and 
S*^ depending on whether the the defect is in its ground 
(g) or excited state (e). In the following, we will take a 
row/column index equal to one (two) for the left (right) 
electrode so that the diagonal (off-diagonal) elements cor- 
respond to reflection (transmission) amplitudes (see Fig- 
ure P). We choose the scattering states to be the eigen 
states of the system when the defect is in its ground state 
and the barrier is totally reflecting hence: 5? = Sij. S'^ 
is an arbitrary unitary matrix with reflection probability 
R = \Sfi\'^ < 1. We will assume that a negative potential 
—V {V > 0) has been applied to the left electrode with 
respect to the right electrode. 

The spectral function, p(lu,V), for absorption by the 
local level is given by 0: 

/oo 
x{tf,V)e''^'fdtf (1) 
OO 

x{tf,V) = {0\U{tf,Om- (2) 

Here |0) is the ground state wavefunction of the complete 
system (the filled Fermi seas in the two electrodes and 
the defect in its ground state), while U{tf, 0) is the time- 
evolution operator for the system between t = and 
t = tf with the defect in its excited state. x(t/,0) is 
the same as the core hole Green's function computed in 

Before discussing the full non-equilibrium case, we 
briefly review the known equilibrium results. When 
V ~ the response of the system is that of the core 
hole problem in a non-separable potential considered in 

miii 

logx(t/,0) = ~i{Eo - A{0))tf ~(3\ogitf^o (3) 

where /3 = 2 (ir) ■ Here e~'^''j are the eigen val- 

ues of S'^. The threshold is shifted from Eq, the en- 
ergy separation in the two-level system, by A(0), which 
is the shift of the ground state energy of the two Fermi 
seas when the scattering defect is in its excited state. 
This standard equilibrium result ^ is well understood 
in terms of the low-lying particle-hole excitations created 
by the rapid switching of the potential, with the princi- 
pal contributions to the logarithm in iPJ from excitations 
with frequencies between and ^o- 




FIG. 2: Typical absorption spectra, p{ijj), in arbitrary units 
for the out-of-equilibrium device sketched in Fig0 The spec- 
tra depend on Sfi = VRe""' and SI2 = Vle*"^ (seelHlandEJ , 
where R is the reflection probability, (a) (01,02) ~ (0,1.5) 
with R = 0.9 (broken line) and R = 0.8 (solid line). As the 
reflection probability decreases the spectrum broadens but 
retains its asymmetric shape, (b) R = 0.5 with (ai,a2) = 
(0,1.5) (solid curve) and (01,02) = (1.5,0) (broken curve). 
The spectrum is sensitive to the difference in ai and 02 . The 
difference between the two curves would show up as differ- 
ences in the spectra on reversing the bias. 

When a voltage is applied across the barrier with the 
defect in its excited state and i? ^ 1, a current will flow 
and the system will become dissipative. For </ <C V~^, 
the spectral response is dominated by excitations with 
frequencies uj V, involving states which do not sense 
the potential drop across the barrier. As a result xi^f: ^) 
is unchanged from its value in equilibrium. 

When tf ^ V^^, the response is controlled by elec- 
trons within the band of width V about the mean Fermi 
energy. We find that 

\ogxitf, V) = -i{Eo - A{V))tf - /3'log (Vtf) + D (4) 

Here the function A(t^) is given by: 

AiV)^r riog(^n(i^)),^ 
J-00 27ri Jq 2m 

(5) 

This expression ((Sj for the (in general complex) energy 
shift of the two Fermi seas, when the defect is in its ex- 
cited state, can be thought of as the generalization of 
Fumi's theorem 0, 0| to the out-of-equilibrium case. 
The exponent /?' in (0J is given by 

P'^Y. (-log(^l,■)/27^^)^ (6) 

The constant term D gives the contribution from excita- 
tions with frequencies between V and which do not 
sense the potential drop across the barrier. To logarith- 
mic accuracy |l6j : 

D=f3\og^o/V. (7) 

Writing S^j = /Re*"^ and comparing the forms for 
P and /?' in Q and ((HJ, we see that the quantity 
— log {Sjj)/2i = —aj/2 + i(logi?)/47r is acting as a com- 
plex phase shift. Its real part, —aj/2, characterizes the 
scattering in the j'th electrode and in (0J describes the 
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effect of particle-fiole excitations in the band of width 
V from the Fermi energy. Its imaginary part (log R) / An 
relates to the lifetime of the excitation. 

The absorption spectrum is found from the Fourier 
transform of x(iy,y) in (Q. Measuring to from luq = 
Eo - Rc(A(y)), it is given by [13 : 



_e-/320n sin (/3> - {p[ - - P'2 log n). 



(8) 

Here we have defined flexpicl)^ = uj/V — z(logi?)/47r 
and written f3' = P[ + if]!^. While the dependence on (3[ 
reflects the total overall scattering on the two sides of 
the barrier as in equilibrium, f3'2 is proportional to the 
difference in the phases of the two reflection amplitudes 
Sfi and 5*12 and its appearance in |(HJ is entirely an out- 
of-equilibrium effect. 

When R = 1, the term multiplying in ((SJ is pro- 

portional to the theta function 6{iu) and describes the 
usual sharp threshold in p{u!). With i? < 1 it leads to 
a smearing of the threshold (see Figure |2J). As pointed 
out in 1^, this broadening of the threshold reflects the 
existence of 'negative energy excitations' in the system 
involving a hole in the left electrode and a particle in 
the right electrode. From an experimental point of view, 
the below threshold broadening with its functional de- 
pendence on the phases of the reflection amplitudes and 
its overall energy scale fixed by the bias are probably the 
key signatures of the non-equilibrium effects we are de- 
scribing. The sensitivity to the difference in scattering 
phase shifts (this difference is proportional to p!^) would 
show up in changes in the line shape on reversing the bias 
and should also be observable. 

The derivation of the overlap x(t/) follows quite closely 
that of MA j^. We introduce the operators ai{e) which 
annihilate particles on the i'th side of the barrier with 
energy e in eigen states of the system with the defect in 
its ground state {S = 1). The effect of the time-evolution 
operator U acting between t = and tf on states aj|), 
where |) is the true vacuum with no particles, is given by 



Ual{e)\)^Y. I de'<J,,{e,e')a]{e')\). (9) 



One can show that for states near the Fermi energy (see 
[Tsf for example) a is given by: 

1 

a,, (e, e') = e"'^"*/ — / S,, {t)e'^'-' ^'dt (10) 
provided that the adiabaticity condition 

is satisfied. In © S{t) = 5^ for t < and t > tf and 
S{t) = S'^ for Q < t < tf and we have suppressed the ex- 
plicit dependence of S on energy. When computing the 



low frequency asymptotics, this becomes a slow depen- 
dence on (e -t- e')/2, and can be neglected. 
The overlap x{tf) can be written 



x(t/) = (0|t/|0)==dctV 



(12) 



where the prime indicates that the operator determinant 
is to be taken only over the occupied states in the two 
filled Fermi seas. This reduces in the equilibrium case 
to the determinant in 0. With zero chemical potential 
in the right electrode and treating the (non-equilibrium) 
Fermi distribution as the diagonal operator /ij(e, e') = 
e')d{—{e + V{2 — i)) allows us to write 



Xitf) = det(l-/ + /a) 
logx(</) = Tr(log(l-/ + /a) 
= C{V,tf) + Ti-f\oga 



/logc^) 



(13) 
Tr/loga 
(14) 



where the operator determinant is now the full determi- 
nant taken over all states and the trace, Tr, is the trace 
over energy and channels. The last term in the expres- 
sion (|14|l can be found by explicitly carrying out the in- 
tegral in pO|l . This gives that (Ty(e, e') = SijS{e — e') — 
Xij{e — e'). The logarithm can then be expanded as 
a power series in the matrix X After evaluat- 

ing X" term by term and then resumming we obtain: 
Tr/log cr = -i{Eo - A(0))i/ + {Vtf /2TTi) (log S)n. The 
difference between this and —i{Eo — A{V))tf in ||SJ) is 
contained in the function C(tf, V). 

To evaluate C{V,tf) we introduce S{t, A) where 



S{t,X) =exp{\\ogS{t)), 



(15) 



so that S{t, 1) = S{t). We now apply the following gauge 
transformation: 



^a(6) 
S{t, A) 



a(e,t) = 
Sit, A) 



e^^^*a(e) (16) 
= e*^^*5(i,A)e-'^^* (17) 



Here L is the diagonal matrix with Ln = 1 and L22 — 0. 
This has the advantage of eliminating the chemical poten- 
tial difference between the two electrodes at the expense 
of an added time-dependence for S when t S [0, t/]. After 
switching to the time-representation (in which the trace, 
Tr, becomes a trace over channels and an integral over 
time) and substituting for a from l(Tn|) . C{tf,V) can be 
written 



C{tf, V) = Tr dX ((1 - / + fSr'f - fS'') ^ 



(18) 



Using a parallel argument to that of we find that 



(1 - / + fS)-' = Y+ ((1 - /)y-i + JYI 



(19) 



where Y± = Y{t ± zO, A). Here Y{z,X) is an analytic 
(matrix) function of complex z in the complement of the 
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cut along the real axis between z = and z = tf, and 
satisfies: 

y_y+ ^ = S{t,X) and Y{z,X) const for \z\ ^ cx3. 

(20) 

If there is no tunneling between electrodes {S^ diagonal), 
this matrix RH problem can be shown to be the same 
as the homogeneous part of that solved in 0. After 
substituting H19(l into H18|l . using the fact that in the 
time-representation (after the gauge transformation [TH|) 
f{t, t') = i(27r(t-i'+z0))"i and letting t' ^tio compute 
the trace, Tr, we finally obtain 



1 /•*/ 

d\ / 

"'0 



\ dt + 



-1 ^-1^ 



dt. 
(21) 

Here tr denotes a trace over channel indices. 

Solving for x(tj,y) is equivalent to solving for the 
quantity Y{z,X). For small V, we can expand the ex- 
ponential factors in S{z, A) fsee ll7|l as e**^^ = 1 ± iVz. 
In this case 



Y[z, A) = exp 



1 

27ri 



log 



log5(z,A) 



(22) 



solves the RH problem. For \z\ — > oo, the exponent (and 
hence Y) tends to a constant as required. If T^i/ <C 
1 we can insert this result into (|21|l and compute the 
integrals over t and A. This yields the equilibrium result 
Although there are corrections to the equilibrium 
iy = 0) solution for Yj^ which are linear in Vt^ these 
cancel out after taking the trace in (|21|l . Corrections to 
C{tf,V) can therefore only be of order {Vtf Y' or higher. 

For times t/ > a general solution to this type of 

matrix RH problem is not known. The form H22() for Y+ 
is still valid for < t < V"^ and tf > t > tf - V^^. The 



integral over times close to the branch points of Y then 
gives the contribution varying as D = log(fo/^) in O- 
However, although the form for Y in H22() still satisfies 
the discontinuity condition along the cut, the exponent 
is unboimdcd for large \z\ and hence (|22|l is useless as a 
starting point for solving for y+ for t ^ V~^. Following 
the derivation of Hi, we find that: 



Y+{t,X) 



ip+{t, A) when t <0 

J •^)^ ^,^(t, A) when 0<t<tf 

?/'+(i, A) when tf<t 

(23) 

is asymptotically correct for t ':s> V ^. Here 7(i, A) = 
5i2(t,A)/5ii(i, A) and i^+it.X) =^{t + iQ,X) where 



■0(2, A) = exp log ■ 



log5n/52*2 



'-To 



log (511^2*2) 



-7"3 



(24) 

The corresponding function Y{z^ A) is not analytic across 
vertical cuts in the complex z-plane through the points 
z = and z — tf, with discontinuities which decay as 
e^^l^l or e~^l^~*J'l. (These factors show that we can- 
not describe the reverse bias case by taking F < in 
(|23|l . Instead Y+ takes a different form for negative V .) 
After inserting the solution H23|) in (|21|l and computing 
the integrals over A and t, we obtain the first two terms 
in Q). The term obtained after differentiating 7 in H23() 
and adding to the term from Tr/ log a in (|14|) , leads after 
some algebra to the term —i{Ea — /S.{V))tf. Differenti- 
ating {tf,X) in H23|) leads to the term proportional to 
log Vt f . The constant term is derived using the form H22() 
for valid for small t and t — tf as discussed above. 
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